Chiral confinement in quasirelativistic Bose-Einstein condensates 
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In the presence of a laser-induced spin-orbit coupling an interacting ultra cold spinor Bose-Einstein 
condensate may acquire a quasi-relativistic character described by a non-linear Dirac-like equation. 
We show that as a result of the spin-orbit coupling and the non-linearity the condensate may become 
self-trapped, resembling the so-called chiral confinement, previously studied in the context of the 
massive Thirring model. We first consider ID geometries where the self-confined condensates present 
an intriguing sinusoidal dependence on the inter-particle interactions. We further show that multi- 
dimensional chiral-confinement is also possible under appropriate feasible laser arrangements, and 
discuss the properties of 2D and 3D condensates, which differ significantly from the ID case. 
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Although cold gases are typically neutral, artificial 
electromagnetism may be induced by several means, in- 
cluding rotation [T] , manipulation of atoms in optical lat- 
tices [2113111], and the use of laser arrangements [5|. Inter- 
estingly, seminal experiments on optically created gauge 
fields have been recently reported [Bj. Artificial electro- 
magnetism has attracted a growing attention in recent 
years, partially due to the possibility of achieving non- 
Abelian gauge fields [H [5], which establish fascinating 
links between cold gases and high-energy physics |H|7l[8]. 
A striking example is given by the possibility of inducing 
quasi-relativistic physics in cold atoms despite the ex- 
tremely low velocities involved 0. In particular, under 
proper conditions cold atoms may experience an effec- 
tive spin-orbit coupling, which leads to a Dirac cone in 
the dispersion [TD], resembling the case of yet another 
paradigm of modern physics, namely graphene [TTl I12j. 
Similar phenomena are expected in cold atoms and 
graphene including Veselago lensing [3 [131 [S] • 

Interparticle interactions lead to inherent nonlinear- 
ities in Bose-Einstein condensates (BECs). At suffi- 
ciently low temperatures the BEC physics is described 
by a nonlinear Schodinger equation similar to that found 
in nonlinear optics [L,. Resemblances between both 
fields have been successfully explored in recent years, 
most remarkably in what concerns the physics of soli- 
tons [13 [ini dZl [13 [T3] , for which nonlinearity and dis- 
persion compensate leading to a non-dispersing solution. 
Nonlinearity plays also an important role in high-energy 
physics. Indeed, non- linear Dirac equations (NLDEs), 
and more generally non-linear spinor fields, have been 
studied extensively, starting with the pioneering works 
of Ivanenko [20 , Weyl [21 , and Heisenberg [22,. These 
equations may present also localized solutions [23l |24| . 

In this Letter we explore the nonlinear physics of a 
multicomponent BEC (also called spinor BEC pS ) in 
the presence of an optically-induced spin-orbit coupling. 
In the low-momentum limit, the spinor BEC may be de- 
scribed by a particular type of two-component NLDE. We 
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FIG. I: One possible laser arrangement for creating a spin- 
orbit coupled gas, based on an atomic tripod configuration. 

show that for the case of attractive interactions, the con- 
densate may become self-trapped, resembling chiral con- 
finement, previously discussed in the context of the mas- 
sive Thirring and Gross-Neveu models [33l|34|. Contrary 
to two-component vector solitons, discussed in nonlinear 
optics and cold gases [33 |3B] , the two components are not 
only coupled by the interatomic interactions, but by the 
spin-orbit coupling as well. This leads to remarkable new 
features of these self-localized BECs, which present in 
ID a peculiar sinusoidal dependence with the interaction 
strength. Furthermore, whereas with usual short-range 
interacting BECs solitons are stable only in ID, we show 
that 2D and even 3D chiral confinement becomes pos- 
sible by means of feasible laser arrangements, although 
the properties of the self-trapped BECs largely depend 
on the dimensionality of the system. 

In the following we consider a Bose gas of atoms with 
an accessible internal tripod scheme (Fig. 1) formed by 
three ground F = 1 states {mp = 0,±1) and an ex- 
cited F — mp — state (this may be the case e.g. of 
the transition 55i/2(F = 1) ^ 5P3/2(F = 0) in ^^Rb 
or the transition 2'^ 6*1 ^ 2'^Po in '*He*). The ground- 
state levels with mp = —1, and 1 are linked to the 
excited state by means of three lasers with, respectively, 
CT4., TT, and cr_ polarizations, Rabi frequencies f2i^2,3(r), 
and phases 6*1,2,3 (r)- The atom-hght interaction results 
in the appearance of two dark states |-Di.2(^), which are 
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linear superpositions of the three lowest bare states. If 
n = (X^i^D^^^ is large compared to any other energy 
scale, including two-photon detuning, Doppler and Zee- 
man shifts and interaction energy, we can neglect transi- 
tions out of the dark state manifold. Assuming that the 
atoms are loaded into this manifold we may express any 
general state as |^'(r,i)) = ^.^^^ ^'^(r, t)|£',i(r)), where 
^'i(r,t) is the wave function for atoms in \Di). By using 
the Schrodinger equation for the tripod scheme, includ- 
ing the atomic motion, and projecting onto the dark state 
manifold, one finds the effective Schrodinger equation [5 
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where 'i>^ = (^'i,^'2)- Here p denotes the mo- 
mentum operator and m is the atomic mass. Note 
that the spatial dependence of the laser arrangement 
leads to an effective vector potential matrix (also 
called the Mead-Berry connection [571 [35]) A„m = 
ih{Dn{r)\'\/Djn{r)). In addition, V„„ = Vn{r)Sn,m and 
= (SV2m)(7?„(r)|VB(r))(S(r)|Vi^™(r)) are effec- 
tive scalar potential matrices, where l-B(f^) is the so- 
called bright state, i.e. the linear combination of the 
ground levels which does couple with the excited state. 
Note that Vn can include any detuning of the n-th laser 
from the resonant transition. 

The exact form of A„„, $„m and Vn depends on the 
incident laser beams and their relative phase and inten- 
sity ratios. To obtain non-trivial gauge potentials the 
light fields need to be shaped carefully. There are many 
techniques to shape light fields in 2D, such as using 
shaped apertures or algorithms to calculate the phase- 
hologram patterns which shape the beam in a chosen 
transversal plane |29j . Shaping arbitrary light beams in 
3D is more challenging but certainly possible [M]. For 
3D shaping holographic techniques are often used [31] . 
This results in a remarkable flexibility to achieve dif- 
ferent gauge fields (for more details see Ref. [5]). In 
this paper we are particularly interested in externally 
inducing an effective spin-orbit coupling. This may be 
achieved by employing three co-propagating lasers along 
e.g. the z axis (5*1 — S2 = S3 — kz), with constant 
I rial and spatially dependent transversal profiles fii = 
jrial cos ^(cc, ?/)e*'^^, = \^3\sm(j){x,y)e'^''^ . Assuming 
a density modulation along x such that 0(r) = \/2kx, 
we obtain the effective spin-orbit coupling A = hnayCx, 
which we employ in our discussion of the ID localized 
solutions. If the density modulation has however a po- 
lar symmetry on the xy plane, i.e. 0(r) = %/2kp (with 

= x'^ + y"^), then A = HnayCp, which we consider for 
the case of 2D chiral confinement. Finally, setting the 
laser detunings A such that Vi = V2 = hA — H^k^ /2m, 
and V3 = —Vi — 2hA, one obtains for both ID and 2D 
arrangements €> + V = HAa^ ■ 

The spin-orbit coupling term p • cr, where cr = 
^y£x,p leads to a single-atom dispersion law charac- 



terized by two branches E±{p) = {p^ + /i^k^)/2to ± 
1/2 

{p?A'^ + H'^k'^pI p/m?^ . The i?_ branch presents two 
minima (or a continuous ring of minima) . As a result the 
ground-state of the many-body system may become frag- 
mented, as recently discussed in Ref. (3^. Fragmentation 
(and the consequent absence of condensation) would pre- 
clude the use of the Gross-Pit aevskii (GP) formalism em- 
ployed below for the analysis of ground state properties. 
The analysis is however well justified if we consider an 
initial (well defined) scalar BEG in one of the ground 
state levels, and adiabatically switch on the lasers, which 
allows a transfer of the BEG (in absence of dissipation) 
into the dark-state manifold. 

Inter-atomic interactions do play an important role in 
the properties of the Bose gas at low temperatures. In the 
following we consider the case in which the interaction 
energy ^ Ml, and hence we can still restrict ourselves 
to the dark-state manifold. Note that the ground states 
\j = —1,0,1) constitute a spin-1 Bose gas. Short-range 
interactions are dominantly s-wave and due to symmetry 
may occur only in two different channels with a total spin 
of the pair and 2, which are characterized by the corre- 
sponding s-wave scattering lengths oq and 02 |25j . These 
scattering lengths are in principle different, although in 
practice they are very similar. Below, and for simplicity 
of the discussion, we consider — a2 — a (for oq ^ 02 
the system remains in the dark-state manifold, but the 
equations coupling the two dark-state components are by 
far much more complicated). In the following we intro- 
duce the coupling constant g ~ 4,TTlT?a/m. 

Within the GP formalism, the interacting bosons in 
the dark-state manifold are described by a spinor GP-like 
equation with a spin-orbit coupling, which in its time- 
independent form acquires the form 
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(p + hua-y + hAa^ + g^^ ■ ^ 
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where /i is the chemical potential, and we consider g < Q. 
For a wavepacket with {p) = and momentum width 
Ap <C 2hK, \j2hAm we can safely neglect the p^ term in 
Eq. ([2]). This results in a nonlinear Dirac-like equation. 



p • cr -I- (7^ - 7^*^ • ^ # 
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where e — /i/SA, 7 — jgl/SA^g for a dimension d, and we 
have employed as length unit = hn/mA. As discussed 
below, the solutions which follow have a typical width 
> Zo- As a consequence, the neglection of the p^ term 
above is valid only for sufficiently small detunings A ^ 
2/iK^/m. 

As already mentioned NLDEs have being extensively 
studied in high-energy physics. Interestingly, it has been 
shown, in the context of chiral confinement in the massive 
Thirring and Gross-Neveu models, that Eq. ^ supports 
in the ID case an exact self- localized solution [331 [51]. 
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In the following we show that these ID self-trapped solu- 
tions present remarkable properties, which have not been 
explored in the high-energy context due to the limita- 
tions of the original physical models. In addition, we 
shall show that multi-dimensional chiral confinement is 
also possible with feasible laser arrangements, although 
it differs significantly from the ID case. 

We study first a ID scenario along the x-direction, as- 
suming a sufficiently strong transversal harmonic con- 
finement (with frequency u)±) in the yz plane, such that 
hijj± ^ fj, (we assume that the temperature is also much 
lower than the transversal oscillator energy). We con- 
sider that the system experiences a spin-orbit coupling 
(see above) p ■ a = PxO'y- Using ^ = ri{cos(p, sinip)'^ 
Eq. (|3| transforms into two coupled equations: 

dip 

dx 
drj 

dx 



— e — cos(2(^) + ^rf 
= — ?7sin(2(^). 



(4) 
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which lead to rf [e^ + ^rf /2 — cos2ip) = const. Im- 
posing localization (77 ^ for x — *■ ±00) one obtains 
£ + ^rj^ /2 — cos2(p — 0. Inserting this conservation law 
back in Q and ([s]) we obtain the localized solution [55] 
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where we have introduced the notation (3 — (1 — e)/(l+e) 
and A = — e'^. From these expressions it is clear that 
only solutions with jej < 1 (i.e. \E\ < HA) are localized. 
Imposing the normalization of the ID wavefunction 
dxrf = 1 leads to the condition 7 — 4tan~^ yT?- As 
a consequence both normalization and localization con- 
ditions are only fulfilled for < 7 < 2tt. In this regime 
the energy presents a remarkable sinusoidal dependence 
£ = cos(7/2). Note that contrary to the high-energy 
case, where the constraint to positive energy fermion 
states demands 7 < tt |33j . in ultracold gases 7 > tt 
is also possible. This leads to a peculiar behavior of the 
BEC wavefunction when 7 approaches 27r. In the regime 
< 7 < 27r the square width of the soliton is 
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12sin'*(7/2) 
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which clearly diverges both at 7 = and 2-k. The di- 
vergence has however a rather different meaning in both 
cases. For 7 = there is no localized solution, as one 
could expect due to the absence of interactions. On the 
contrary for 7 = 27r, the localized solution exists but 
acquires a Lorentzian form 1/(1 + 4a;^), with divergent 
(x^) but finite half-width at half-maximum 1/2. Direct 
numerical simulations of Eq. (1) confirm the existence of 
these localized solutions, which remain well self-localized 
for hA < 0.1hK'^/2m. 
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FIG. 2: Energy e as a function of the interaction strength 7 
for a 2D and 3D arrangement. Note the existence of a solution 
for all 7 > 4.36 in 2D and a double solution for all 7 > 11.94 
in 3D. 

Due to the form of ip{x), the relative distribution of 
the population between the two components also presents 
a peculiar dependence in 7. The chirality x (defined 
as the difference between the populations in \Di) and 
\D2)) becomes X = / dxrj^ cos{2ip) = |sin(7/2), which 
monotonously decreases from 1 (all atoms in \Di)) at 
7 = to (equal admixture of both dark states) at 
7 = 27r. The density profile of both components is also 
quite different, since \D2) presents a node at a; = 0, 
whereas \Di) is maximal at the center. 

Note that the discussed sinusoidal dependence has its 
origin in the ID normalization of (|7]). We should hence 
expect a rather different behavior in higher dimensions. 
This is indeed the case. We consider in the following a 
2D scenario on the xy plane (we assume a strong con- 
finement along z) in the presence of a spin-orbit coupling 
p • cr = puy, where = Px+Py- Eqs. (6]) and ^ can be 
readily generalized to the 2D case by replacing x by the 
radial coordinate p. The normalization of ^ is however 
different 27r J prf'{p)dp = 1, which requires 
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sin 2s 



(9) 



where e — cos(2s) (0 < s < 7r/2) and L(s) = 
— /p ln(cos s')ds' is the Lobachevsky function [37]. By in- 
verting ([9]) we obtain £(7) (Fig.[2|. As mentioned above a 
localized solution requires |e| < 1, which is only possible 
for 7 > 7c = 27rln2 ~ 4.36. Therefore, contrary to the 
ID localization which may occur for < 7 < 27r, the 2D 
localized BEC requires a minimal interaction strength 7c 
at which the BEC width diverges. In addition both local- 
ization and normalization conditions are simultaneously 
fulfilled for arbitrary 7 > 7c- For increasing 7 > 7c, e de- 
creases monotonously from 1 to —1. For 7 > 15, e ~ — 1 
and the localized wavefunction converges to a Lorentzian 
shape Tf{p) = 87r7-i(l + 4p2)-i. Fig. [s] shows 2D lo- 
calized BECs. Note that the behavior of the chirality x 
in 2D differs from the ID case. In particular for e ^ 1, 
X ^ 1 as in ID, and for e —1, x ^ —l, i-s--, I-D2) 
dominates. Note also that as in ID the \D2) state has a 
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FIG. 3: Left column: e = -0.99 and 7 : 
total density rf is shown at the top. 



= 18.63. The localized 
The density of each 



component ^l/i 2 are strikingly different, with component 1 
forming a peak in the center. Most particles are occupying 
component 2. Right column: e = 0.99 and 7 = 4.36. The 
component 1 is forming a peak in the centre where now most 
particles are located. The x and y coordinates are in units of 
Iq and the densities in units of l/^o- 

minimum at the center, whereas \Di) is there maximal. 

Finally, a similar solution may occur also in 3D if p-u = 
pcTy, where = p'^+Py+pl- In that case the 3D normal- 
ization requires that 7 = 2s(7r^-4s^)/3sin^(2s) (Fig. [2]). 
The localized solution exists only for 7 > 7c ~ 11.94. 
Contrary to 2D where the BEC width diverges at 7c and 
hence the wavefunction experiences a smooth crossover 
from localization into delocalization, in 3D at 7 = 7c, 
£ < 1 and the BEC width remains finite. Hence in 3D 
there is an abrupt transition between localization and 
delocalization regimes. Finally, let us point out that, in- 
terestingly, for 7 > 7c, there are actually two localized 
solutions. At 7 cxD, one of the solutions becomes un- 
bound and the other a Lorentzian shaped function. 

In summary, the interplay between interactions and 
an optically induced spin-orbit coupling may lead to self- 
localized Bose-Einstein condensates in ID, but also in 2D 
and 3D. Self- localization in NLDE was previously studied 
in the context of chiral confinement in massive Thirring 
and Gross-Neveu models. However, for the case of cold 
gases a novel regime of parameters (and dimensionalities) 
are possible, allowing for remarkably rich physics which 
largely depends on the system dimensionality, ranging 
from a sinusoidal interaction dependence in ID, to two 
possible self-trapped solutions in 3D scenarios. These re- 
sults provide exciting new perspectives for the nonlinear 
physics of condensates in artificially-induced gauge fields. 
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